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O ■ Abstract 
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q^ | Motivated by the strong astronomical evidences supporting that huge black- 

holes might inhabit the center of many active galaxies, we have studied the 
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integrability of oblique orbits of test particles around the exact superposition 



of a black-hole and a thin disk. We have considered the relativistic and 



the Newtonian limits. Exhaustive numerical analyses were performed, and 

bounded zones of chaotic behavior were found for both limits. An intrinsic 

relativistic gravitational effect is detected: the chaoticity of trajectories that 

do not cross the disk. 
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In recent years, strong observational evidences have supported that huge black-holes, with 
masses between 1O 6 M and 1O 1O M , might inhabit the center of many active galaxies [|TJ. 
These evidences have motivated many investigations on the black-hole-disk system. Some 
approximate |2j and numerical || results were obtained, and in Qj some exact axisymmetric 
solutions describing systems containing the superposition of non-rotating black-holes and 
static thin disks are presented and discussed. Such systems have no net angular momentum, 
and a possible explanation for their stability is that the disk particles move under the action 
of their own gravitational field in such a way that there are as many particles moving to 
one side as to the other || . This counterrotating interpretation has been frequently used to 
describe true rotational effects (see, for spherical systems, and, for cylindrical systems, 

0)- 

Here, we consider the integrability of oblique orbits of test particles in the exact static 

black-hole-disk system. There are examples in the literature of chaotic motion involving 
black- holes: in the fixed two centers problem || , in a black-hole surrounded by gravitational 
waves U , and in several core-shell models with relevance to the description of galaxies (see 
for a recent review). We mention also that geodesic motions in some static axisymmetric 



spacetimes were considered in [JTTJ] , where chaotic behavior was detected for several solutions 
containing N-point Curzon-like singularities. 

We are manly interested in bounded motions close to the black-hole, so one assumes 
that the disk is infinite and has an associated homogeneous Newtonian density. We have 
also considered the Newtonian limit of the black- hole-disk system, i.e., we have also stud- 
ied the integrability of oblique orbits of test particles under the Newtonian gravitational 
potential corresponding to the superposition of a monopole and an infinite homogeneous 
disk. Bounded zones of chaotic behavior were found for both the relativistic system and 
for its Newtonian limit. We also notice that the chaotic regions in the relativistic cases 
are typically larger than in the corresponding Newtonian ones. This is a consequence of an 
intrinsic relativistic gravitational effect we found: the chaoticity of trajectories that do not 
cross the disk. 



Let us start by the Newtonian limit of the black-hole-disk system. The equations of 
motion for test particles in this case are very simple. We use cylindrical coordinates (r, 6, z) 
with the monopole, with mass M, located in the origin. The disk corresponds to the plane 
z — 0. The gravitational potential in this case is given by 

M 
V(T,B,z) = —== + a\z\, (1) 

where a is a positive parameter standing for the mass density of the disk. The angular 
momentum L in the z direction is conserved, and we can easily reduce the three-dimensional 
original problem to a two-dimensional one in the coordinates (r, z) and with the Hamiltonian 
given by 

TT f 2 z 2 M L 2 

H = 1 ; H -+az, (2) 

2 2 Vr 2 + z 2 2r 2 ' K J 

for z > 0. For z < 0, H is obtained by substituting a — * —a. The Hamiltonian (0) is 
smooth everywhere except on the plane z — 0. Moreover, the parts of the trajectories 
restricted to the region z > (or z < 0) are integrable. The Hamilton- Jacobi equations for 
a two-dimensional system with the potential 

„, , M L 2 

v ^ = -W77^w +az (3) 



can be separated in parabolic coordinates |12| , leading to the second constant of motion 

^ n r2 •/ • -\ a i Mz r „z , AS 

C = R z - a— = r(zr - rz) r z -\ . + L — , (4) 

2 V J 2 Vr 2 + z 2 r 2 ' V ; 

where R z is the z component of the Runge-Lenz vector 

M 
R= - T ===(rr + zz) + v x L, (5) 

Vr + z z 

where L stands for the total angular momentum. In this case, with the two constant of 
motions H and C, the equations for the trajectories of test particles can be reduced to 
quadratures in parabolic coordinates [0. We can take advantage of these results to study 
the trajectories of the original system (g). 



Suppose that a bounded trajectory (H < 0) starts somewhere in the region z > 0. The 
time evolution of the trajectory will be governed by the integrable potential (^) until the 
plane z = is reached, say at r = ri, f = r\, and z — Z\ < 0. After crossing the disk, the 
trajectory will be also governed by a potential like (^), but now with a negative parameter 
a. The trajectory evolves until the plane z = is reached again, at r = r 2 , r = r 2 , and 
£ = z-i > 0, when the signal of a must be changed again in order to get the correct time 
evolution, and so on. The relevant Poincare's section in this case will be given by the 
intersections {(ri,ri), {rz,fz), . . .}. Since the motion in the regions z > and z < are 
integrable, any irregularity in the oblique trajectories shall be credited to the signal changes 
of a between the integrable regions. In particular, trajectories that do not cross the disk are 
integrable, in contrast to the relativistic case, as we will see latter. 

The trajectories in the integrable regions can, in principle, be described by means of 
Jacobi elliptic functions. The expressions are considerably complicated, and we have aban- 
doned the hope of constructing analytically the Poincare's sections. Nevertheless, we could 
construct them very accurately by solving numerically the system in the integrable regions 
and matching the trajectories appropriately on the disk plane. In the Figures 1 and 2, we 
show typical Poincare's sections across the plane z = 0. Fig. 1 presents a low-energy situ- 
ation (H = — 0.4, L = M — l,a — 0.1) where the integrability seems to be preserved in a 
large region, while Fig. 2 presents a section (H = — 0.15, L — M — 1, a — 0.1) revealing a 
widespread chaotic behavior. We could obtain thousands of intersections for each trajectory 
with a cumulative error, measured by the constant H, inferior to 10~ 12 . 

A key point is that the Poincare's sections showed in Fig. 1 and 2 are present for any 
finite values of M and a. In particular, we will have wide zones of chaotic motion, as 
presented in Fig. 2, for any finite values of M and a. This is a consequence of the invariance 
of the equations of motion under the transformations 

r -> Xr, M -> M, z -> Xz 
a -► X- 2 a, t -► X Z 'H , L -► \^XL, (6) 



H-> \- l H 

and 

r -> AV, M -> A'M, 2 -> A'z, 
a^A ,_1 a, *->A't , L^X'L, (7) 

A > and A' > 0. Poincare's sections are invariant, up to constant rescalings, under such 
kind of transformation, and hence it is possible to choose suitable values for L and H in 
order to have, for instance, the Poincare's section of Fig. 2 for any finite values of a and M. 
Now, let us focus the relativistic static black-hole-disk system. We start with the Weyl 
metric [[HJ describing static axisymmetric spacetimes, 

ds 2 = - e *lr,z) d1 2 + e u{r,z)-v{r,z) ^2 + ^ + ^ e ~u(r,z)^ (g) 

The metric describing the superposition of a black-hole with mass M and a thin disk with 
homogeneous associated Newtonian density a can be obtained from the metric of a black- 
hole immersed in a pseudo-uniform field, which is described by the potentials 

R 1 + R 2 - 2M 



u(r, z) = az + In 



R 1 + R 2 + 2M'' 



uj{r,z) = -—r + In ^-^ + a{R 2 -R 1 ), (9) 



where R\ = J(M — z) 2 + r 2 and R 2 = J{M + z) 2 + r 2 . The first terms in the right-handed 
side of (H) correspond to the pseudo-uniform gravitational potentials, which can be obtained 
from the Curzon metric by a limiting process [ 14|| . They are called pseudo-uniform because 



their Newtonian limit corresponds to a uniform gravitational field. The potentials (|9]) satisfy 
Einstein vacuum equations everywhere except at the origin r = z = 0. 

One can transform our metric in Weyl coordinates (t, r, 6, z) into spherical ones (t, r§, 0, <p) 
by doing 



r§ — 2Mr s sin <p, z = (r s — M) costp. (10) 



For a = 0, this coordinate transformation puts the line element (g), with the potentials 
(0), into the Schwarzschild form. From ([T0|), one sees directly that the Schwarzschild radius 
r$ = 2M corresponds in Weyl coordinates to the rod r = 0, —M < z < M, and that 
the circular equatorial photonic orbit at rg = 3M, (p = ir/2, corresponds to the orbit at 
r = V3M, z = 0. We note also that under the transformation (|T0|), any infmitesimally thin 
disk in Weyl coordinates remains infmitesimally thin in spherical coordinates. For more 
details on the geometry of the Weyl metric, see [H| and the references therein. 



Our black-hole-disk system corresponds to the juxtaposition of two black-hole-pseudo- 
uniform field solutions: the solution (^) for z > and a solution with a reversed uniform field 
(a — ► —a) for z < 0. The obtained metric has C° components and obeys Einstein vacuum 
equations everywhere except in the plane z = 0. The only nonvanishing components of the 
energy-momentum tensor are 

u{r,z)-u){r,z) / M \ 

u(r,z)-ui(r,z) yr 

»» = T -' = " 4, VM* + / M - (11) 

Our solution obeys the weak energy condition (e > 0) everywhere. The counterrotating 
velocity V of the particles of the disk is given by V 2 = pee/e |[J]. For a relativistically 
consistent solution, one usually demands that V 2 < 1. For the present case 

V= , M . (12) 



We have V^ 2 < 1 for r > y3M. If we change from the original Weyl to Schwarzschild 
coordinates, we find that it corresponds to ro, > 3M. Hence, we conclude that our solu- 
tion is indeed relativistically consistent, since it is a well-known result that there are no 
circular orbits inside the photonic orbit (rs = 3M). We are obtaining the consistent result 
that the centrifugal-gravity balance inside the photonic radius can only be maintained for 
superluminal velocities [|J. 

Now, we need the geodesic equations for the black-hole-disk system. Since one always has 
two independent Killing vectors for static axisymmetric spacetimes, the geodesic equations 

6 



for the Weyl metric (§) can be cast as 



r + /i(r, z) (f 2 - i 2 ) + 2/ 2 (r, z)ri + #i(r, z) = 0, 
^ + /2(r, 2) (i 2 - f 2 ) + 2/i (r, z)ri + # 2 (r, 2) = 0, 

where the dots denote derivation with respect to s. For the potentials ©, we have: 



(13) 



fi{r,z) =r 



1/R 1 + l/i? 2 1/Rj + l/Rl 



a. 



R 1 + R 2 + 1M 



1 

R~ 2 

ff . I 1/Ri + I/R2 l/R\ + l/Rl a ( 1 1 
H\ r i z ) — z 



2 V-Ri 



T 



M 



R X + R 2 + 2M 2 

r 1/R 1 -1/R 2 1/Rj-l/R 2 . 



2 V i?x i?2 
a / 1 1 
2 \R\ R 2 



h 2 (r,z) 



gi(r,z) = hi(r,z) 
g 2 (r,z) = hx(r, z)h 2 (r, z) 
where 



R l + R 2 + 2M 2 

8M(R 1 + R 2 )r AR±R 2 L 2 



(i?i + R 2 + 2M) 2 r 3 

SM(R x {z + M) + j^z - M)) 
(i?i + R 2 + 2M) 2 



exp(2az) 



+ 2a.R1.R2 



a 

~ 2 

R x + R 2 - 2M 
{R 1 + R 2 + 2Mf 

R x + R 2 - 2M 



(14) 



R Y + R 2 + 2M 



a 



hi(r, z) = exp —r + a(Ri - R 2 ) 



h 2 {r,z) 



E 2 



+ 



exp(2az) 



(15) 



(i?i + R 2 - 2M) 2 r 2 (R x + R 2 + 2M) 2 
The expressions (JTJ]) and (|I5|) are valid for z > 0. For z < 0, /1, / 2 , g>i, and (72 are obtained 



by doing a — >• —a. The constants £/ = —e u( - r,z H and L = r 2 e~^ < - TV " l # are, respectively, 
the constants of motion associated to the Killing vectors d/dt (the energy) and d/dg (the z- 
angular momentum). Besides of these two constants, the requirement of a timelike trajectory 
leads to a further one 



Qab% •& 



-1 



^(r,z)-u(r,z) U2 + j,2\ + f^Ov) _ -gl^v^z) 



(16) 



We stress here that the system (|T^)-(|T6|) corresponds to the simplest case of superposition 
of black-holes and disks. For the superposition with a non- homogeneous or a finite-radius 
disk, the metric tensor is considerably involved ||, and even the task of writing down the 
geodesic equations for oblique orbits seems to be unrealistic. 



We could solve flT3|) numerically by using the same scheme used in the Newtonian case. 
However, we point out a crucial difference: the parts of the trajectories contained in the 
region z > (or z < 0), in contrast to the non-relativistic case, are themselves non integrable. 
In other words, the trajectories of test particles around a black-hole immersed in a pseudo- 
uniform gravitational field are non integrable. This case should be added to the list of 
known examples of systems presenting qualitatively distinct behavior in their relativistic and 
Newtonian limits |ll|, stressing the deep differences between both theories. This fact clearly 



indicates that the relativistic black-hole-disk system is more chaotic than its Newtonian 
limit, and we could indeed verify that the chaotic regions in the relativistic limit are typically 
larger than in the corresponding Newtonian one. 

Figure 3 presents a Poincare's section revealing a widespread chaotic behavior for the 
trajectories around a black-hole immersed in a pseudo-uniform field (E = 0.975, L = 3.8, 
M = 1, and a = 5 x 10~ 4 ). We could obtain thousands of intersections for each trajectory 
with a cumulative error, measured by the constant (|T6|), inferior to 10~ 12 . Since the rela- 
tivistic equations are invariant under only one type of rescaling, our conclusions concerning 
the universality of the Poincare's section presented in the Fig. 3 are more limited than in 
the Newtonian case. It is not a surprise to find a strong chaotic behavior of the oblique 
orbits around the black-hole-disk system, since they are obtained by matching properly on 
the disk plane the (chaotic) trajectories of two black-hole-pseudo-uniform field systems. 

We finish discussing our implicit assumption that there is no interaction between the test 
particle and the disk. It was assumed that the test particle may cross the disk many times 
without suffering any trajectory deviation. This can be an acceptable hypothesis if one 
considers the disk particles as very small when compared to the test particle and that we do 
not have many crosses. However, even in the case of very small disk particles, some trajectory 
deviation is expected to occur after a huge number of intersections. We have considered the 
following case of a weak interaction between the disk and the test particle: each time the 
test particle crossed the disk, its velocities (r, z) and L were perturbed by small and random 
amounts such that H (or E in the relativistic case) was maintained constant. Qualitative 

8 



equivalent Poincare's sections were obtained. However, a more realistic interaction should 
also permit small changes in H (or in E). In other words, it would be worth to study the 
behavior of our models under more general stochastic perturbations ||16[ . These points are 
now under investigation. 
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FIGURES 




FIG. 1. Poincare's section (r,r) across the plane z = for oblique orbits, with H = —0.4 and 
L = 1, around the superposition of a monopole with mass M = 1 and an infinity homogeneous 
disk with surface density a = 0.1. 
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FIG. 2. Poincare's section (r, r) across the plane z = for oblique orbits, with H = —0.15 and 
L = 1, around the superposition of a monopole with mass M = 1 and an infinity homogeneous 
disk with surface density a = 0.1. 
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FIG. 3. Poincare's section (r, r) across the plane z = for orbits with I? = 0.975 and L = 3.8 
around the exact superposition of a black-hole with mass M = 1 and a weak pseudo-uniform 
gravitational field (a = 5 x 10 -4 ). 
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